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We consider dipolar interactions between heteronuclear molecules in low-dimensional geometries. 
The setup consists of two one-dimensional tubes. We study the stability of possible few-body 
complexes in the regime of repulsive intratube interaction, where the binding arises from intertube 
attraction. The stable dimers, trimers, and tetramers are found and we discuss their properties 
for both bosonic and fermionic molecules. To observe these complexes we propose an optical non- 
destructive detection scheme that enables in-situ observation of the creation and dissociation of the 
few-body complexes. A detailed description of the expected signal of such measurements is given 
using the numerically calculated wave functions of the bound states. We also discuss implications 
on the many-body physics of dipolar systems in tubular geometries, as well as experimental issues 
related to the external harmonic confinement along the tube and the prospect of applying an in-tube 
optical lattice to increase the effective dipole strength. 

PACS numbers: 67.85.-d,68.65.-k,42.50.-p 



I. INTRODUCTION 

Ultracold polar molecules with anisotropic long-range 
interactions have generated a lot of interest recently. In 
recent experiments heteronuclear molecules could be pre- 
pared in their rotational and vibrational ground-state 
and cooled to temperatures close to quantum degener- 
acy [1-7 , where many exotic many-body states have 
been predicted [HJ [9] . In a three-dimensional sample 
one faces the problem that attractive head-to-tail in- 
teractions of the dipoles can destabilize the system by 
strong particle loss [10] . In order to overcome this prob- 
lem, molecules can either be dressed by ac-external fields 
leading to strongly repulsive interactions at small inter- 
particle distances or optical lattices such that particles 
cannot approach each other in the head-to-tail direc- 
tion [TTHT4] . Recently, experiments have entered the 
quasi-two-dimensional regime by applying an optical lat- 
tice potential to the three-dimensional sample and ob- 
served large geometric effects on the loss rates [15] . While 
this system does not provide the idealized stack of two- 
dimensional traps, it spurs hope that this can be realized 
in the near future. In the meantime a number of theo- 
retical works on layered dipolar systems have appeared. 
For a single layer of dipolar molecules, p-wave superfluids 
P~6j [17] and density-waves [T8H2T] have been predicted. 
For large dipole moments and perpendicular dipoles a 
single layer is expected to crystallize in a triangular ar- 
ray [8, 22-28 analogous to the famous Wigner crystal 
in the Coulomb gas [29j [30] . For bilayers or more gener- 
ally multilayers, the long-range character of dipolar forces 
causes interlayer interactions that can lead to interlayer 



coherence and ferromagnet-type states [31], BCS-BEC 
crossover phenomena [32, 33 for bilayers, and dimer- 
ized pairing [34] for multilayers. Interlayer interactions 
have been studied extensively for analogue systems with 
long-range Coulomb interaction in quantum Hall [35] or 
graphene bilayers [36 and in electron-hole semiconductor 
bilayers where exciton superfluids or supersolids might 
be realized [37j [38] . The few-body physics in single- or 
bilayer systems have also been studied recently [39H45] , 
and some aspects of the many-body problem in coupled 
one-dimensional geometries have been discussed [46H5T] . 
The details of Efimov physics for three dipolar bosons in 
three dimensions have also been discussed [52 . In a re- 
cent letter [53] we predict the stability of complex bound 
states of dipolar molecules, that consist of more than one 
molecules per tube. In the present work we provide the 
details of our results and elaborately analyse the full pa- 
rameter space wherein complexes are stable. 

Our basic setup consist of two parallel one-dimensional 
tubes at a fixed distance containing ultracold polar 
molecules. While the motion is only along the tubes, 
the long-range character of the dipole-dipole interaction 
means that the molecules interact across the tubes. We 
consider which few-body bound states for the bi-tube 
system are stable as the direction and the strength of 
the dipoles are varied. We focus on the regime where 
intratube interactions are repulsive so that the binding 
interactions of the few-body states comes from the in- 
tertube attractions. The formation of a trimer in our 
setup strongly resembles the formation of charged exci- 
tons in semiconducting nanotubes [54] which are also ID 
structures. A positively charged exciton consists of two 
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FIG. 1: (Color online) (a)- (d) Schematic pictures of the few- 
body complexes labeled by the number of molecules in the 
tubes A and B. (a) 1-1 dimer (b) 1-2 trimer, (c) 2-2 four- 
body state and (d) 1-3 four-body state. Also indicated are 
the coordinates of the individual particles. 



holes that are bound together by an additional conduc- 
tion electron. Analogously, in our setup two molecules 
in the same tube repel each other, but can be bound 
together by the attraction from a molecule in the other 
tube. 

We label the few-body complexes by the numbers Na- 
Nb of particles in tubes A and B that participate in the 
complex as shown in Fig{l] Below we discuss the stability 
of the 1-1, 1-2, 2-2 and 1-3 complexes in detail. Then we 
propose an optical measurement scheme to detect vari- 
ous complexes for which a schematic picture of the corre- 
sponding experimental setup is shown in Fig. [2] We show 
that the light scattering is strikingly different for all con- 
sidered few-body states. Since it is non-destructive it 
allows the in-situ observation of the creation and dissoci- 
ation of the few-body complexes. We also discuss various 
extensions of our work including the case of parallel layers 
instead of tubes, multi-layer or multi-tube systems and 
effects of few-body bound states in the many body prob- 
lem corresponding to a finite particle density. Technical 
details of the calculations are given in the Appendix. 



II. MODEL 

The general Hamiltonian for N dipoles is given by 
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where z, j label the particle and run from to N — 1. We 
assume that all dipoles are aligned by an external electric 
field and describe strength and direction of the dipoles 
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FIG. 2: (Color online) Visualization of setup. Molecules move 
in two parallel tubes named A and B that are located in xy- 
plane and are separated by distance A. Red circles illustrate 
dipolar molecules. The dipole moment is represented by an 
arrow d = d(cos(#) cos(</?), cos(#) sin(<^), sm(#)). <p denotes 
the angle between tube and dipole within xy-plane and ft the 
angle out of the plane. The optical detection scheme consists 
of an additional probe beam and a detector and is discussed 
in Sec|V] 



by the vector 

d = d(cos(^) cos((£>), cos($) sin((^), sin(#)). (2) 

A schematic picture of the corresponding experimental 
setup is shown in Fig. [2] The interaction between two 
dipolar molecules is given by 
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with D 2 = cP/Aireo and cos (p r d = r • d/ (rd), where r de- 
notes the relative position ri — T2 of the two molecules. 
The motion of the molecules is reduced to two ID tubes 
at positions (x,y,z) = (x, ±A/2,0) where A is the in- 
tertube distance. We consider deep ID optical lattices 
where intertube tunneling can be negelected and where 
the transverse confinement length, is much smaller 
than the distance A between the tubes. In the limit 
A ^> l± the intertube interaction, Vi(x), and the in- 
tratube interaction, Vq(x), depend only on the distance 
between the particles along the tubes in the x-direction. 
However, we take into account that the intratube inter- 
action is modified at small distances by the transverse 
part of the wave function [26] 

D 2 

V (x) = -— (1 -3cosVcos 2 ^)A 3 /o(Ax), where (4) 
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with x = x/A, and where l± is the transverse con- 
finement length. In the strict ID limit one replaces 
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FIG. 3: (Color online) (a) Main panel shows the intratube 
potential Vo(x) for perpendicular dipoles ip — 90°, $ = 0° for 
various strengths of the transverse confinement A see Eq. Q. 
Tilting the dipoles leads to a suppression of Vb without chang- 
ing its functional form. As shown in the inset the intratube 
interaction vanishes at the magic angle ipM- (b) The inter- 
tube interaction for ip = ipM and ip = 90° both with $ = 0°. 
Since particles in different tubes are always separated by at 
least the intertube distance A the dependence of Vi(x) on 
the transverse confinement is of order 1/A 2 and is neglected 
in this paper. Energy in units of h 2 /mA 2 and length in units 
of A. 
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with x = x/A. Figure [3] shows the intratube and in- 
tertube interaction and illustrates their dependence on 
tilting angle and transverse confinement. 

The characteristic energies for dipole interaction and 
kinetic energy are D 2 / A 3 and h 2 /mA 2 respectively 
which defines a dimensionless measure for the interaction 
strength given by 



U = mD 2 /AH 2 



(7) 



Uo determines the competition of interaction to kinetic 
energy. The classical limit where the kinetic energy of the 
bound state can be neglected corresponds to Uo — » 00. 
If not stated, otherwise we measure length in units of A 
and energies in units of h 2 /mA 2 . 

We consider few-particle complexes of up to four 
molecules in total as shown pictorially in Fig. [I] We char- 
acterize the configurations by the number of molecules in 
tubes A and B. For example the four body complex 1-3 
consists of 1 molecule in tube A and 3 in tube B. The 
interaction potentials for the various complexes are given 
by: 
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X 3 fo(Xx) — >> 1/x 3 in Eq. Q. The finite transverse con- 
finement reduces the intratube repulsion at small dis- 
tances (of order l±) and regularizes the 1/x 3 behav- 
ior in the strict ID case [26]. Since Vb(0) = (1 — 
3 cos 2 (p cos 2 1?) ^- and the radius of the core is of 

order Zj_, the short range behaviour of the intratube in- 
teraction is similar to a contact interaction 5(x)D 2 
The short range part is particularly important for bosonic 
molecules. For A ^> 1 this short-range repulsion exceeds 
intertube attraction and bosons become hard-core and 
behave as fermions in ID. However, this is different at 
tilting angles where the prefactor (1—3 cos 2 (p cos 2 #) = 0. 
For $ = 0, this condition defines the so-called magic an- 
gle, ip M = arccos(^). 

Corrections to the intertube interaction in quasi- ID 
scale with 1/A 2 and as we have used A = 5 and 10 in 
our calculations these can be neglected. In turn, the 



with Since the interaction only depends 

on the relative coordinates, their center of mass mo- 
tion can be separated. To a given set of parameters 
(/?, A, [To, Na, Nb we determine the eigenspectrum by 
solving the Hamiltonian of the relative motion numer- 
ically using a finite difference method which turns the 
solution of the Schrodinger equation into a diagonaliza- 
tion problem of large sparse matrices as discussed in Ap- 
pendix[Aj The stability of a few -body complex is checked 
in two complementary ways. First we confirm that the 
eigenenergies of the N = Na+Nb particle state is smaller 
than any other smaller few-body state, so that there is 
an energy penalty against dissociation. For example a 
stable 1-2 state has to have smaller energy than the 1- 
1 state. This argument is strictly true only in the zero 
density limit and in absence of an external parabolic con- 
finement potential, so that it is always possible to have 
a free particle with zero energy. The second method is 
to check the average distance between the particles. In a 
bound state all distances have to be finite. 
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FIG. 4: (Color online) Different regimes for few-body states 
in tubes. In the left upper corner above the dashed-dotted 
line f dxVi(x) < and there is always a bound dimer state. 
In the center of the figure between dashed-dotted and dashed 
line there is intertube attraction but J dxV\{x) > 0, so that a 
dimer only exists above a critical interaction strength. To the 
right of the (red) dashed line V\ (x) > everywhere. In the left 
lower corner limited by the solid (black) line, the intratube 
interaction is repulsive. 



III. FEW-BODY BOUND STATES 

In this section we work in the low density limit and 
discuss how few-body states consisting of a small number 
of particles in each tube depend on the system parame- 
ters Uq, (/?, $ and A as well as on whether the molecules 
are fermions or bosons. Effects related to a finite den- 
sity per tube will be addressed in Sec. 



IV We focus 



on the regime where the intratube interaction is repul- 
sive so that few-body complexes are bound due to the 
intertube attraction. In the case of intratube attraction 
higher transverse bands or even chemical reactions might 
become important and the effective intratube interaction 
given in Eq. Q should be modified accordingly at small 
length scales. 

The angles $ and cp introduced in Fig (2] will determine 
the ratio between intratube and intertube interaction. 
For a given (p the intertube attraction V\ is strongest 
for i? = and in order to stabilize larger complexes we 
will mainly focus on $ = in this section. An exemp- 
tion is the dimer which turns out to be always stable at 
d = 0, but can be made unstable at d > as shown 
in Fig. [4j The intratube interaction Vq becomes repul- 
sive for cp > arccosf 7= 1 , rr ). At 1? = the intratube 

r v V3cos('iy) / 

interaction becomes repulsive for cp > cpM- 



A. Dimer 

We start our discussion with the simplest few-body 
state, the dimer, that consists of one particle per tube 
as sketched in Fig. [IJa). The stability of the dimer as 



function of the angles $ and cp is shown in Fig. [4j For 
t? = 0a dimer can always profit from the attractive part 
of the intertube interaction and is stable for all cp and 
Uq. This can be shown by noting that J V\-i(x)dx < 0, 
which is sufficient to ensure a bound state in ID. Even 
in two dimensional systems where J V\-i(x)dx = a 
dimer always exists [42 . However, in ID the dimer 
can be made unstable by rotating the dipoles out of the 
x — y plane corresponding to d > 0. For general d the 
integral over the intertube interaction is J dxVi(x) = 
2 (cos 2 ((p) cos 2 (tf) - cos(2tf)) D 2 /A 3 . There are three 
regimes: 1) For j dxV\{x) < there is a dimer bound 
state for any Uq. This regime is in the upper left cor- 
ner and is limited from below by dashed dotted line in 

which corresponds to (p 



Fig. 



For JdxVi(x) > but Vi(x m i n ) < there is a bound 
dimer above a critical interaction strength. Here x m i n is 
the position of the minimum of V\x around x = (as seen 
for instance in Fig. [3^b)). This regime is in the center 
of Fig. [4] and is limited to the left by the dashed dotted 
line and to the right by the dashed line. The dashed line 
corresponds to d = arccos(l/\/3). In the third region 
V\(x) > for all x and there is never a bound dimer. 
This situation occurs to the right of the dashed line in 
Fig. [4j We note that for a single dimer it does not mat- 
ter whether the molecules are fermions or bosons as long 
as tunneling between the layers can be neglected. Fi- 
nally, the intertube interaction is only weakly dependent 
on the transverse confinement strength for A ^> 1 as the 
correction to the strict ID limit scales as A -2 . 



B. Trimer 

Next we discuss the trimer, which consists of one par- 
ticle in one tube and two particles in the other tube as 
depicted in Fig. [TJb). As is the case of the dimer, the in- 
tertube attraction is responsible for binding the trimer. 
For a trimer to be stable, the energy gain from inter- 
tube attraction has to overcome the energy cost from 
intratube repulsion and from localization. In this sub- 
section we consider d = where the intertube attrac- 
tion is strongest. The stability of the trimer can be in- 
creased by tilting the dipoles close to the magic angle 
LpM = arccos(l/v // 3) or by increasing L7 . The latter re- 
duces the cost of localization energy and the former re- 
duces the intratube repulsion. Figure [5|a) shows how the 
trimer becomes stable as the interaction strength Uq is 
increased at a tilting angle cp = 57° close to the magic an- 
gle. Above a critical interaction strength of Uq = 2.5 the 
energy of the trimer is smaller than three free molecules 
as well as one free molecule and one dimer. Experimen- 
tally Uq can be tuned by the strength of the homogeneous 
electric field. However, in order to reach large values of 
Uq, which are needed here, one needs to choose heteronu- 
clear molecules with a large electric dipole e.g. 6 Li 133 Cs. 

The instability of the trimer at small interaction 




FIG. 5: (Color online) (a) Energy of the dimer state 1-1 (solid 
(black) line) and the trimer state 1-2 (dashed (blue) line) for 
(p = 57° and A = 10. Vertical gridline at Uo = 2.5 indicates 
critical interaction strength above which trimer is stable. For 
the parameters shown here bosonic and fermionic trimer have 
almost the same energy and larger few-body complexes are 
unstable, (b) Spatial extension of dimer and trimer as defined 
in Eq. (12). Energy in units of h 2 /mA 2 length in units of A. 
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FIG. 6: (Color online) (a) Contourplot of the the trimer po- 
tential (|9| for cp = 57° as function of the relative coordinates 
ui — (—2xo + xi + X2) / V§ and U2 — (xi — X2) / V% in units of 
A, see Eq. (A4). Intralayer repulsion gives a large potential 



barrier at U2 — 0. Brighter areas correspond to higher values, 
(b), (c) Controuplot of the wave function for fermionic and 
bosonic trimers. The bosonic wave function is almost given by 
magnitude of the fermionic wave function which is expected 
for hard core bosons. Here Uo — 10 while the remaining pa- 
rameters are as in Fig. [5] 



strength Uo is also visible in the spatial extension of the 
trimer. Denoting by x n the position of particle n and 
by X c = TV -1 x n the center-of-mass, where N is the 
total particle number, we define the spatial extension of 
a complex as 

z2 = ^E(^- x c) 2 = ^£(^-* m ) 2 . (12) 

n m<n 

Fig. [5jb) shows the spatial extension of the dimer and 
trimer. The dissociation of the trimer into a free particle 
and a dimer at small Uo leads to a diverging the spatial 
extension of the trimer for Uo < 2.5. We note that the 
energy of the trimer smoothly approaches the energy of 
the dimer. 

For the tilting angle <p = 57° and the interaction 
strengths shown in Fig. [5j the energy of the trimer is al- 
most independent of whether the molecules are bosons or 
fermions. The reason is the strong intratube repulsion at 
short distances, which turn bosons into hard core bosons 
that behave like fermions in ID. This can be appreciated 
in Fig. [6] The intratube repulsion causes a large poten- 
tial barrier at U2 = as shown in Fig.[6ja). As shown in 
Fig.[6^b) and (c), the wave function of the bosonic trimer 
is to a good approaximation given by the magnitude of 



the wave function of the fermionic trimer. 

At the magic angle, the effective intratube interac- 
tion vanishes and the difference between trimers made 
up of either bosonic or fermionic molecules is maximal. 
Fig. [TJa) shows the potential, Vi_2, of the trimer and 
the wave functions for bosonic or fermionic molecules. 
The relative wave function for the fermionic and bosonic 
trimer are shown in Fig. [7^b) and (c). For fermions, 
the antisymmetry of the wave function under particle ex- 
change enforces that the wave function vanishes at = 
although the potential is minimal there. Due to the an- 
tisymmetry of the wave function the fermionic trimer is 
only stable above a critical interaction strength even at 
the magic angle. In contrast, the wave function of the 
bosonic trimer has a maximum at U2 = for ip = <pm 
as shown in Fig. [7^c) and the bosonic trimer forms at 
arbitrarily small L^o at the magic angle. 

Fig. [8] shows the stability of the trimer states both for 
fermionic and bosonic molecules as function of the tilt- 
ing angle ip for transverse confinement strengths A = 5 
and 10. The solid black line shows the energy of the 
dimer which has little dependence on the tilting angle 
for the regime shown. The thick dashed line shows the 
energy of the fermionic trimer at a strong transverse 
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FIG. 7: (Color online) Same as Fig. [6] but for cp = ipM, where 
intratube interaction vanishes, (a) Contourplot of the trimer 
potential shows a minimum at U2 = 0. (b) Fermionic trimers 
has a node at ui = due to Fermi- statistics, (c) Bosonic 
trimer has maximum at U2 = 0. 



ergy is minimal since intratube repulsion vanishes. By 
increasing tp, the intratube repulsion increases and at 
about ip = 60.5° the trimer becomes unstable. The 
thick dash-dotted line shows the energy of the fermionic 
trimer at a smaller transverse confinement A = 5. Low- 
ering the transverse confinement reduces the intratube 
repulsion at short length scales and thus lowers the en- 
ergy of the trimer. At the magic angle, the intratube 
vanishes exactly and within our model (which assumes 
A ^> 1) the energy of the trimer is independent of A. 
For the fermionic trimer, the transverse confinement has 
a rather small effect, since the short length repulsion is 
less important due to the Fermi-statistics that forbids 
two fermions to be at the same point. However, as we al- 
ready discussed in Fig. [7] bosonic molecules can share the 
same single-particle wave function and therefore bosonic 
trimers have a much lower energy than fermionic ones for 
small intratube repulsion. The energy of bosonic trimers 
is shown by the thin lines in Fig. [8] The energy of the 
bosonic trimers increases much faster than for fermionic 
trimers as cp is moved away from the magic angle. Addi- 
tionally, the transverse confinement has a stronger effect 
on the energy for bosonic molecules than for fermionic 
molecules. However, away from the magic tilting angle 
the intratube repulsion is sufficiently strong to suppress 
the overlap between the two molecules in the same tube 
and the bosonic and fermionic trimers have nearly the 
same energy. This situation is shown in Fig. [6j 
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FIG. 8: (Color online) Energy of dimer and trimer few-body 
states as function of the tilting angle ip > ipM for fixed inter- 
action strength Uo = 10. The subscripts F and B indicate 
fermionic or bosonic molecules. For strong lateral confinement 
(A ^> 1), the dimer energy is independent of particle statis- 
tics and on the lateral confinement. In contrast, the trimer 
energy depends on both A and statistics, in particular close 
to the magic angle ipM ~ 54.73° where intratube interaction 
is strongly suppressed. At the magic angle only particles in 
different tubes interact with each other and this interaction is 
almost independent of A. Note that we compare total energy 
of dimer and trimer although they consist of different number 
of particles. Energy in units of h 2 /mA 2 . 



confinement of A = 10. At the magic angle its en- 



C. Four-body states 

There are two possible four-body states, the 2-2 state 
with two particles in each tube and the 1-3 state with 
one particle in one tube and three particles in the other 
tube. They are shown schematically in Fig.[TJc) and (d). 
Similar to the trimer, the stability of these states is lim- 
ited to tilting angles close to the magic angle <pm and 
sufficiently large interaction strength Uo. Again this is 
needed so that the attractive intertube interaction over- 
comes the repulsive intratube interaction and the cost 
of localization. For the parameters shown in Fig. [5j the 
four-body states are unstable. However, closer to the 
magic angle and for sufficiently strong interaction the 
four-body states are stable. 

Fig. [9] shows the stability of the 1-3 states both for 
fermionic and bosonic molecules as function of the tilting 
angle cp and for different transverse confinement strengths 
A. Fig. [9ja) shows the fermionic case. The thick lines 
correspond to a transverse confinement A = 10. For 
(p < 56.3° the 1-3 state (thick dashed line) is stable and 
has a lower energy than the 1-2 trimer (thick solid line). 
For cp > 56.3° the 1-3 complex dissociates in a trimer 
and a free particle. The dissociation can also be seen 
in the spatial extension of the 1-3 complex (not shown 
here) which diverges for cp > 56.3° in a similar way as 
for the trimer in Fig. [5] The thin lines in Fig. [9^a) corre- 
spond to A = 5 and show qualitatively the same behav- 
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(b) U = 10 Stability of bosonic 1 -3 
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FIG. 10: (Color online) Energy of four-body states in the 
2-2 configuration and the energy of two independent dimers 
as function of the tilting angle ip > ipM for fixed interaction 
strength Uo — 10. While the dimer energy is independent 
of A and statistics, the 2-2 state depends on both. We note 
that this time we compare complexes with the same number 
of particles. Energy in units of h 2 /mA 2 . 
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FIG. 9: (Color online) Energy of trimer and four-body states 
in the 1-3 configuration as function of the tilting angle ip > 
ipM for fixed interaction strength Uo — 10. Both 1-2 and 1-3 
states depend on lateral confinement A and statistics. The 
stability of the 1-3 state is limited to a small interval of tilt- 
ing angles close to the magic angle ipM ~ 54.73°. Again we 
note that we compare total energy of 1-2 and 1-3 complexes. 
Energy in units of h 2 /mA 2 . 



ior. Fig. [9[b) shows the bosonic case. We note that for 
bosons the complexes have a much lower energy than for 
fermions at the magic angle (fM- However, the energy 
rapidly increases for larger tilting angles. For bosons, 
the short distance behavior of the intratube interaction 
is relevant close to the critical angle and there is a strong 
dependence on the transverse confinement. 



Fig. 10 shows the stability of the 2-2 states both for 
fermionic and bosonic molecules as function of the tilting 
angle cp and for different transverse confinement strengths 
A. At the magic angle the 2-2 state has lower energy than 
two independent dimers. Increasing intratube repulsion 
by increasing the tilting angle the 2-2 states soon become 
unstable. Notice again that bosonic states have a much 
smaller binding energy than fermionic ones. 



D. Stability diagram 

The regions of stability of the various complexes in the 
Uo-ip plane is shown in Fig. 11 Above the (blue) dashed 
line only the dimer is stable, while below it the trimer 



also becomes stable. The (green) dotted line marks the 
position at which the 1-3 state becomes stable and be- 
low the (red) dashed-dotted line the 2-2 state is stable. 
The diagram shows that with increasing number of par- 
ticles in the complex the stability regime is reduced to 
a small intervall starting at the magic angle. Further- 
more, by comparing the stability diagram for fermionic 
and bosonic molecules a major difference is evident. For 
fermions, there is a critical interaction strength even in 
the absence of intratube repulsion, while bosons do not 
have a critical value. 



Our results on the the stability of the few-body states 
at the critical angle are collected in Fig.[l2]and[l3j As we 
can see complexes with two or more fermionic molecules 
in the same tube are only stable above a critical inter- 
action strength even in the absence of intratube interac- 
tion. Due to Fermi-statistics the wave function has to 
have a node for x\ = x^ which increases the kinetic en- 
ergy and increases the localization energy. In contrast, if 
the complex is made up of bosonic molecules, the wave 
functions can be finite at x\ = x^ and there is no critical 
interaction strength. The regimes of stable trimers is ex- 
tended for both bosons and fermions and lies in the range 



<Pm < ip < 62° for the values of Uo shown in Figs. 12 and 
[l3l This is a rather broad interval and implies that the 
study of trimers should be relatively easy. However, for 
the tetramer states the window of angles is much smaller 
and extends only at most one or two degrees away from 
LpM- Complexes larger than trimers therefore need a 
higher degree of fine-tuning of the setup and are more 
difficult to access. 
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(a) Fermion stability diagram 
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(b) Boson stability diagram 
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FIG. 12: (Color online) Energy of few-body states 1-1, 1-2 and 
1-3 in units of J^ A2 at the magical angle (fM as function of 
Uo for A = 10. (a) Even in the absence of intratube repulsion 
the 1-2 and 1-3 configurations of fermonic molecules are only 
stable above critical interaction strength indicated by vertical 
lines, (b) For bosonic molcules, however, the 1-2 and 1-3 
configurations are stable at any interaction strength at ip = 



FIG. 11: (Color online) Stability of different complexes as 
function of interaction strength Uo and tilting angle cp for 
A = 10. Each region is labeled by its stable complexes. For 
large tilting angle, only the dimer is stable but as one ap- 
proaches the magic angle ipM, other complexes become sta- 
ble. Below the (blue) dashed line the trimer is stable, below 
the (green) dotted line the 1-3 state is stable and below the 
(red) dashed-dotted line the 2-2 state is stable, (a) Fermionic 
molecules. All few-body states except the dimer have a criti- 
cal interaction strength, (b) Bosonic molecules. Close to the 
magic angle the complexes become stable for any interaction 
strength. 



E. Larger Complexes 

For even larger complexes we could in principle con- 
tinue the methods used here, calculate binding energies 
and make the relevant comparisons. However, already 
from the data presented so far we can infer the influence 
of such states. Consider the case of N = 5 for which a 
2-3 or a 1-4 configuration is possible. According to our 
arguments above, the stability regions for larger com- 
plexes in the bitube case will shrink to a small sliver just 



around ip = (fM- In contrast, the three- and four-body 
complexes are relevant over a broader range of parameter 
space, making them more accessible in experiments. 



F. 2D Layers 

When the tubes are replaced by 2D planes the extra 
non-confining dimension increases the cost of localization 
and bound complexes are less stable. However, we note 
that the classically preferred configurations are the same 
for bilayer and bitubes. Therefore trimers or even larger 
complexes will also be stable in bilayers, however, only 
for larger dipole strength [45 j. An interesting difference 
between ID and 2D concerns the dimer. The appearance 
of bound dimers in 2D is guaranteed at arbitrarily small 
couplings as long as the interlayer interaction integrates 
in space to zero or a negative value [43]. In 2D this pre- 
condition is fullfilled for any direction of the dipol es for 
the interaction In contrast, we discussed in Sec III A 
that in ID the interlayer interaction can be made purely 
repulsive by increasing $ > 54.7°. For purely repulsive 
interaction the dimer is of course unstable. 
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(b) = Stability bosons 




FIG. 13: 



(Color online) Energy of few-body states 1-1, 

h 2 



1-2, 



2-2 in units of mA2 at the magical angle ifM as function of 
Uo for A = 10. (a) Even in the absence of intratube repulsion 
the 2-2 configurations of fermonic molecules is only stable 
above critical interaction strength indicated by vertical lines, 
(b) For bosonic molcules, however, the 2-2 configurations is 
stable at any interaction strength at ip = ifM> Furthermore, 
the 1-2 configuration has even slightly lower energy than two 
independent 1-1 configurations. 



IV. MANY-BODY PHYSICS 

At finite density one needs to minimize the free energy 
per particle. The stability of complexes is determined 
by comparing energies per complexes that can consist of 
different numbers of particles, e.g. trimer and dimer ener- 
gies. If the interactions between complexes are neglected 
and the system is at equilibrium at a finite tempera- 
ture, then stable complexes will be populated according 
to their energy per particle. At low temperatures and 
for balanced systems the occupation of trimers will be 
low, since the energy per particle is higher than in the 
dimer. However, for an imbalanced system the formation 
of larger complexes can be favorable even at low temper- 
atures. For instance, if the number of particles in one 
tube is twice that of the other, then the trimer would be 
favored over the dimer as soon as it is stable. For the 
1-3 complex one would need a 1-3 ratio of the number of 
molecules in each tube and so forth. 

For increasing densities interactions between com- 
plexes can no longer be neglected. An effective interac- 
tion between complexes can be obtained using the bound- 
state wave functions calculated here [33] . Another impor- 
tant aspect is the statistics of the complexes [55]. While 



for bosonic molecules all complexes will behave as bosons, 
the character of complexes of fermionic molecules will 
change with the number of particles. Then complexes 
with an odd number of molecules will have fermionic 
statistics, while even number complexes will behave as 
bosons. 

In the dilute limit, we expect that when the bind- 
ing energy of a particular few-body state consisting of 
fermionic molecules exceeds the Fermi energy then this 
state is well-defined. This situation will be realized in 
the strong-coupling limit where L^o is large. The Fermi 
energy depends on the density whereas Uq depends on 
the molecules, the electric field used to align them, and 
on the optical lattice needed to generate the geometry. 
These two quantities can therefore be tuned indepen- 
dently. From the strong-coupling side, the question then 
arises what the effective degree of freedom is. While 
the dimer and four-body can be considered as bosonic 
degrees of freedom, the trimer is still a fermion. This 
means that systems dominated by one or the other can 
have different quantum proporties. By changing the di- 
rection or strength of the dipoles we expect a crossover 
between fermionic and bosonic behavior, that should be 
observable using some of the techniques already proposed 

[nj [ii [Mi [sg. 



V. OPTICAL DETECTION 

The few-body complexes can be detected using light 
scattering. Recently, several nondestructive (in the sense 
of the quantum non-demolition, QND) schemes for mea- 
suring the properties of the many-body states in ultracold 
gases observing scattered light have been proposed [56]- 
[62] . Among them, the method developed in Refs. [56H59] 
is the most relevant to the present system, as it explic- 
itly uses the sensitivity of light scattering to the rela- 
tive position of the particles forming a complex. This is 
due to the constructive or destructive interference of the 
light waves scattered from the different particles. More- 
over, that method can be directly applied for extended 
periodic structures (many equidistantly spaced tubes or 
layers) and many-body systems, which makes the exper- 
imental realization promising. 

We consider the scattering of the probe light with the 
amplitude given by the Rabi frequency Vt p = d^Eo/h (Eq 
is the probe-light electric field amplitude and do is the in- 
duced dipole moment), cf. Fig. [2] To increase the signal, 
the scattered light can be collected by a cavity, and the 
photons leaking from the cavity are then measured. Al- 
ternatively, the measurement of photons scattered can be 
made in a far-field region without the use of a cavity. 

After several assumptions (the small tube radius, far 
off-resonant light scattering, detection in the far field 
zone), the light scattering has a simple physical inter- 
pretation. The scattered light amplitude is given by the 
sum of the light amplitudes, scattered from each molecule 
(cf. Fig. [5J where several molecules are schematically de- 
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picted within two tubes). Each term has a phase and 
amplitude coefficient depending on the position of the 
molecule as well as on the direction and amplitude of the 
incoming and outgoing light waves: 



C 



=A,B 



J dxhi(x)u*(x, pi)u p (x, p^, (13) 



where the sum is over two tubes A and B, hi(x) is the 
operator of particle linear density. For the free space 
scattering, the value of C corresponds to the Rayleigh 
scattering [63]. Adding a cavity to the setup the scat- 
tering is increased and C = —ig s ^t p /(A a tz) with k being 
the cavity decay rate, g s is the molecule-light coupling 
constant, and A a is the light detuning from the reso- 
nance, cf. Refs. [56-59 . In Eq. (13), u p ^ s (x,pi) are the 



mode functions of probe and scattered light at the tube 
positions Pa,b, which contain the information about the 
propagation directions of probe and scattered light waves 
with respect to the tube direction. For the simplest case 
of two traveling light waves, the product of two mode 
functions takes the well-known form from classical light 
scattering theory: u*(x, pi)u p (x, pi) = exp [i(k p — k s )r^], 
where kp 5S are the probe and scattered light wave vectors, 
and is the molecule position. 



Equation ( 13 ) is valid for any optical geometry and can 



describe the angular distribution of the scattered light. 
However, an important conclusion of Refs. 56 59j was 
that some information about the many-body state can 
be obtained even by a simple measurement of the pho- 
ton number scattered at a single particular angle, which 
is enough for our purpose. Moreover, as it was shown, 
the particularly convenient angle of measurement corre- 
sponds to the direction of a diffraction minimum, rather 
than Bragg angle (diffraction maximum). At the direc- 
tions of diffraction minimum any classical (possibly very 
strong) scattering (leading to a signal scaling )is sup- 
pressed, and the light signal exclusively reflects the quan- 
tum fluctuations of the particles. 

We now fix the optical geometry as follows (cf. Fig. 
1). The incoming probe light is a traveling or stand- 
ing wave propagating at the direction perpendicular to 
the tubes, which gives u p (r) = R(x) exp(ik p y) (for the 
traveling wave) or u p (r) = R(x) cos (k p y) (for the stand- 
ing wave) and includes the transverse probe profile R(x) 
of an effective width W. To perform the measurements 
at the direction of a diffraction minimum, the scat- 
tered light is measured along z direction. For the free 
space detection, or the traveling- wave cavity, this gives 
u s (r) = exp(z/c s z), while for the case of a standing wave 
cavity, u s (r) = cos(k s z). Without loss of generality, we 
can assume u s (r) = 1 at the tube position z = 0. The 
absolute values of the wave vectors are equal to their 
vacuum quantities k PjS = 27r/Aii g ht. 

An important property of such a configuration (illumi- 
nation and detection at directions perpendicular to the 
tubes), is that all atoms within two different tubes scat- 
ter light with a fixed phase difference with respect to 



each other, independently of their longitudinal position 
x within the tube. Here we assume a tight transverse 
confinement such that the motion of atoms in the per- 
pendicular directions is frozen. Due to this fact, the aver- 
aging over the probabilistic position of the complex does 
not involve the light phase. Moreover, averaging over the 
probabilistic relative positions within each complex does 
not involve the dependence on the light phase as well. 
At other directions, both those kinds of phase averaging 
are important and would decrease the optical signal and 
the distinguishability of the complex types. Our scheme 
is very convenient since detection takes place at a single 
angle perpendicular to the tubes only and requires no 
scan over multiple angles. 

The operator of the light amplitude reduces to 



C 



(u p (y A )N A (W) + u p (y B )N B (W)) , (14) 



where Na,b(W) are the operators of the effective par- 
ticle numbers in the tubes A and B within the region 
illuminated by the laser beam, 



AW) = f 

J — ( 



fiA y B(x)R(x)dx. 



(15) 



If the laser profile can be approximated by a constant in 
the interval (-W/2,W/2), the operators N a ,b( w ) ex " 
actly correspond to the atom number operators in two 
tubes within the laser beam. 

The classical condition of the diffraction minimum 
is fulfilled, when the expectation value of the light- 
amplitude operator (14) is zero due to the perfect 



cancelation of the expectation values of two terms in 
Eq. ( [l4| ) (i.e. the total destructive interference be- 
tween the scatterers in two tubes). This is achieved 
for u p (y b) / u p {y a) = —{N A )/{N B )- We introduce the 
atom number ratio a = (N A )/(N B ). For the equal 
mean atom numbers (the few-body complexes 1-1 and 
2-2), the optical geometry should be chosen such that 
u s (i/b) /us(va) = —1, which can be achieved if, e.g., 
the tube spacing is the half of the light wavelength, 
A = Aught /2. F° r the few-body complex 1-2, a = 
1/2, and the diffraction minimum is achieved if the 
light wavelength and tube spacing satisfy the condition 
cos(h p y b) / cos(k p y a) = —1/2. This can be achieved, 
e.g., if the position of the tube A corresponds to the 
antinode of the standing wave cos(k p yA) = 1, while that 
of tube B corresponds to k p ys = 2tt/3 or 47r/3, leading 
to the ratios between the tube spacing and light wave- 
length as A = Aught /3 or 2Ah g ht/3. Similarly, for the 1-3 
complex, that ratio can be A « 0.3Aii g ht or 0.7Aii g ht- All 
those example ratios can be indeed larger, taking into 
account the periodicity of the light wave. 

The expectation value of number of photons scattered 
at the direction of diffraction minimum n$ is then given 
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by 



n$ 



(a\a s ) = \C\ 2 \u p (y A )\ 2 x 



(N A (W) -aN B (W) 



(16) 



where u p (i/a) can be easily chosen as 1. This expression 
manifests that the number of photons scattered in the 
diffraction minimum is proportional to the second mo- 
ment of the "rated" particle number difference between 
two tubes in the laser-illuminated region. In general, the 
photon number at the diffraction minimum is non-zero. 
It directly reflects the particle number fluctuations and 
correlations between the tubes. Thus, the number of pho- 
tons reflects the quantum state of ultracold molecules. 

We now apply Eq. (16) to the few-body states dis- 
cussed in this paper. The calculations of the expec- 



tation values in Eq. (16) simplifies strongly if one as- 



sumes the Gaussian shape of the probe light as R(x) = 
exp (— x 2 /W 2 ) with the beam width W. 

Applying Eq. ( 16 ) to the case of one molecule per tube, 



where the formation of dimers is possible, we obtain the 
photon number as 



n 1 



2a 



e 2w 2 



(17) 



where no = \C\ \u p (i/a)\ (W/L)^/it/2 is an effective 
photon number, which can be scattered from a single 
molecule and takes into account the probability of find- 
ing a molecule in an illuminated fraction W j L of the tube 
of the length L. The brackets denote the average over the 
relative wave function. As discussed in Appendix [A] one 
separates the center of mass and writes the two particle 
wave function as \j/(x p, x\) = ^c(uo)^(ui) were uq,u\ 
are defined in Eq.(A2) then 



( X Q- X1 ) 



-I 



du\e w 2 ip 2 (ui). 



(18) 



Equation (17) describes the dependence of the photon 
number on the effective size of the dimer xq — x\. For a 



tightly bound dimer, the expectation value in Eq. (17) is 



(«o-«ir 



1, 



(19) 



while for two independent molecules it approaches 0. 
Correspondingly, in a diffraction minimum (a = 1), the 
light scattering is completely suppressed for a dimer and 
approaches a constant 2no for two unbound molecules. 
This manifests the fact that for a dimer, the molecules 
appear in the laser beam only in pairs and the light scat- 
tered from the two molecules destructively interfere. In 
contrast, for the free molecules, there is a probability that 
only one molecule will be found in a laser beam, while 
another molecule is out of it. In the latter case, the per- 
fect destructive interference is not possible, which leads 
to the nonzero photon number scattered in the diffraction 
minimum. 




FIG. 14: (Color online) Optical detection of a dimer for 
Uo — 1/2 and <p = 90°. (a) Visualization of two different 
regimes depending on the angle (b) Width of the dimer 
given by Eq[l2] in units of A between the dissociating par- 
ticles, (c) Dimer energy in units of J^ A2 ■ (d) Number of 
photons scattered into diffraction minimum, Eq |17| in unit of 
no: Two different photon numbers (zero and 2no) correspond 
to the dimer and two free particles. The light intensity jump 
corresponds to the dissociation (or creation) of a dimer. 



Figure 14 shows the striking sensitivity of the pho- 
ton number to the presence of a dimer or free molecules. 
When the dimer dissociates during tilting of the d angle, 
the photon number sharply jumps from the zero value 
to a value 2no (Fig. 14d). Figure 14 displays the dimer 
size and energy as well as the photon number. The latter 
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shows a characteristic jump as the dimer becomes unsta- 
ble above a critical angle 




FIG. 15: (Color online) Optical detection of trimer for Uo — 
10 and <p = 57°. (a) Visualization of three different regimes 
depending on the angle (b) Relative distances in unit s of 
A between the particles as defined in the text after Eq. ( 20 ) 



showing the consecutive dissociation of the trimer into a dimer 
and unbound molecules, (c) Energy of the state in units of 
■ (d) Num ber of photons scattered into diffraction min- 
imum, Eq. pOl, in units of no: Three different photon num- 
bers (zero, 1/2710, and 3/2no) correspond to the trimer, dimer 
and a free molecule, and three free particles. The light inten- 
sity jumps correspond to the dissociations (or creations) of a 
dimer and trimer. 



For the trimer case 1-2, the photon number reads 



n 



1 + 2a z + 2a l 



-2a ( e 



2W 2 



_ ( x 2~ x l) 



2a ( e 



2W 2 



.(20) 



As in this case, the population of the tubes is imbal- 
anced, the diffraction minimum corresponds to a = 1/2 



and for the tightly bound trimer, the photon number 
in the diffraction minimum is zero. As previously, vari- 
ous terms correspond to the contribution of various rel- 
ative coordinates between the molecules. For the tightly 
bound complex, all distances are small and the expecta- 
tion values (exp( — (xi — Xj) 2 /2W 2 )) « 1, while for large 
interparticle distances they vanish. 

We are interested in three different regimes: 1) Trimer 
is stable. 2) Trimer is unstable but dimer is stable. 3) 
Trimer and dimer are unstable. In order to distinguish 
these regimes we do not only calculate the extension of 
the trimer Z, but calculate the following quantities (where 
the coordinates are given in Fig[T]) 



I 1 = 



max ^2 
mm\2 



r = 







(Max[(x - xi) 2 , (x - 
(Min[(x - xi) 2 , (x - 
((x 1 -x 2 ) 2 /9). 



(21) 

-* 2 ) 2 ]/9>, (22) 
x 2 ) 2 }/9), (23) 
(24) 



Here Iq is the intratube distance measuring the average 
distance between the two molecules in the same tube. 
There are two intertube distances, Z™*™ , , between 
the single molecule in tube A and the two molecules in 
tube B. l™ zn denotes the smaller and \™ ax the larger 
intertube distance. 

Figure 15 shows the relative sizes, energy and photon 
number for the three-body complexes. At small $ the 
trimer is stable, all length are finite and the intensity of 
the scattered light is approximately zero. At $ = 19° the 
trimer dissociates in the dimer and a free molecule. Then 
the larger of the intertube distances Vf 3 ^ diverges and the 
light intensity jumps up, because the condition of the 
perfect destructive interference is not satisfied (the free 
molecule can be within or outside the probe beam). In- 
terestingly, the energy of the trimer smoothly approaches 
the energy of the dimer and therefore the energy shows 
no indication of the dissociation. In contrast the scat- 
tered light intensity jumps to a new plateau. At even 
larger angles At $ > 51° the dimer also dissociates and 
the system consists of three free molecules. Now all dis- 
tances diverge and the light intensity increases further as 
well, manifesting the increased molecule number fluctu- 
ations in the laser beam waist. As shown in Fig. 15, all 
three possible molecular states can be well distinguished 
by light scattering and correspond to three different light 
intensities: zero for a trimer, no/2 for a dimer and a free 
molecule, and 3no/2 for three free molecules. The com- 
plex dissociation corresponds to the jump in the photon 
number. When $ decreases, the creation of the com- 
plexes is manifested by the stepwise decrease in the light 
intensity, which corresponds to the reduced atom particle 
fluctuations. 

Notice that our detection scheme uses the difference 
between plateaus to detect the few-body complexes. In 
a real experiment, there could be additional background 
noise added to the signal. However, the difference in 
the height of the plateaus stays the same. This relative 
measurement of heights is therefore more sensitive than 
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an absolute one, and can help eliminate problems from 
constant noise in experiments. 

In the many-body context, when the density of 
molecules is small, the results for the number of scattered 
photons Eqs. (17) and po] ) (displayed in Figs. 14 'd) and 



|l5[d)). still hold and should be multiplied by the num- 
ber of molecules. The advantage of using the diffrac- 
tion minumum is that the classical coherent scattering 
which scales with the number of molecules squared is 
suppressed, and we get our signal from the fluctuations 
which scale with the number of molecules. We still re- 
quire a low density in order to allow the molecules to 
dissociate to larger distances. In a high density limit, it 
would not make sense to go beyond the average particle 
distance. When this becomes comparable to the com- 
plexes size, we have effectively a crystal configuration 
and no jumps in the signal. 

The analysis of the four-body complex 1-3 can be car- 
ried out in a similar way. It results in three steps in 
the photon number dependence on which correspond 
to the 1-3 complex, 1-2 trimer and 1 free molecule, 1-1 
dimer and 2 molecules, and four unbound molecules. 

While the photon signal from single molecules can be 
small, we have outlined a number of ways above to en- 
hance the signal. This can be done by increasing the 
number of molecules in the tubes, increasing the probe 
laser intensity, the measurement time, or by using a high- 
quality cavity. Additionally, we note that an array of 
many tubes in one plane would also enhance the signal 
proportionally to the number of tubes due to the per- 
pendicular arrangement of probe and detection. Which 
one of these mechanisms is the better option depends on 
concrete experimental setups. 



VI. EXPERIMENTAL ISSUES 

In order to see the complexes discussed here in an ex- 
periment it is the value of Uo that matters. Our calcu- 
lations indicate that the trimer should be stable from 
Uo > 1.5 at the magic angle ip = cpM (see Fig. 12), 
whereas for a larger angle of ip 



57° it is stable for 
Uo > 2.5 (see Fig. [5|. A minumum dipoles strength is 
therefore required to see complexes beyond the dimer in 
current experiments with fermionic polar molecules. The 
strength parameter can be written in the convenient form 



Uq = 0.15 



(t 



m 



amu 



1 ii 

aT 



D 



1 Debye 



(25) 



where 1 Debye is 3.336 • 10 -30 Cm in SI units. For cur- 
rent experiments using 40 K 87 Rb the maximum dipole 
moment is 0.56 Debye [2H5]. Assuming A = 532 nm 
we obtain Uo = 1-12 which is a rather small value for 
investigating few-body complexes. Working with a dif- 
ferent molecule of larger moment therefore seems appeal- 
ing. For instance, a molecule like 6 Li 133 Cs should be able 
to accomodate much larger dipole moments [6] and by 



tuning the external field a large range of Uo should be 
accessible with that combination of alkali atoms. An- 
other advantage of using a much larger dipole moment is 
that the binding energy increases rapidly with Uo- This 
implies that one does not need to be in the degenerate 
regime and can study the bound complexes even with 
thermal samples. 

An interesting alternative idea is to apply an optical 
lattice along the tubes in order to quench the kinetic 
energy of the molecules. This can be understood as in- 
creasing the effective mass, m*, of the molecules over the 
bare value, m*/m > 1. This results in a new effective 
dipole strength, U* = m*U /m > U '. However, the intro- 
duction of such a lattice changes the dispersion relation 
of the molecules in the tubes. The study of bound states 
in periodic potentials with short-range interactions have 
demonstrated modifications of the spectrum as the lat- 
tice depth is increased [64-66 . Similar effects are likely 
for polar molecules when a periodic lattice potential is in- 
cluded in the bound state calculations. Such an elaborate 
extension is, however, beyond the scope of the present 
work and we include only an estimate of the influence of 
a lattice. 

In order to study the fermionic trimer in the 40 K 87 Rb 
system we need an effective dipole strength of at least 
1.5, or even larger to study the dependence on cp. For 
the sake of argument, let us assume that we are inter- 
ested in obtaining an effective U* > 2Z7o- If we assume 
a tight-binding approximation, the effective mass will 
then have to be m* > 2m. Consider a standard opti- 
cal lattice potential of the form V(x) = sEo sin 2 (27nr/A), 
where Eq = 2h 2 7r 2 /m\ 2 is the recoil energy. Numeri- 
cal calculations [67] show that one needs s = 5 to have 
m* ~ 2m when two-body interaction are ignored. The 
tight-binding approximation compares very well to a full 
numerical diagonalization for s > 5 [67] . 

The effect of interactions on the dispersion in the lat- 
tice is small unless the interaction energy considered is 
close to Eq. Let us for simplicity assume that the in-tube 
lattice has the same wavelength as the transverse lattice, 
A = A/2. We then have E = ~ 4.9; 



i.e. 



h 2 

mA 2 ' 



From Fig. [5] we see that the dimer and trimer binding en- 
ergies for Uo ~ 2.5 is ~ 2 ^^2. This implies that the in- 
teraction energy is smaller than the recoil energy and the 
dispersion relation should be given by the tight-binding 
approximation. Another way to check that the complexes 
will not be severely modified by the presence of the lattice 
for small s, is to consider the bandwidth of the lowest lat- 



16_ 5 3/4 e -2^ 



tice band, Wo- It is given by Wo/E 

For s = 5 we find W /E ~ 0.345 or W = 1.7^. The 
latter is comparable to the binding energy for Uo ~ 2.5 
and we therefore do not expect strong modification of the 
bound state energies and structures. 

Finally, we address the external parabolic confinement 
that is always present in experiments with ultracold 
atoms and molecules. This is most often approximated 
by an harmonic oscillator in the plane or tube which is 
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preferably shallow to approximate the homogeneous case. 
For a trapping angular frequency of cj, then the quantity 
of interest is rj = = m ^ -• If rj is small, we 

expect the outer confinement to be negligible. In the re- 
cent experiments at JILA [15 the in-plane external trap 
has a frequency of 36 Hz, yielding 77 ~ 0.13. The exter- 
nal trap is therefore expected to have small influence on 
the bound states. In the opposite limit of large 77, the 
bound complexes must be calculated with the external 
trap included. The methods used will work in that setup 
also. 

We have suggested an optical detection scheme, which 
promises a less destructive, in situ measurement proce- 
dure up to a physically exciting QND limit. The few- 
body complexes can be also observed by other techniques 
including destructive ones. For example, the ID lattice 
depth can be changed periodically. If the shaking fre- 
quency matches the binding energy of a complex it dis- 
sociates thereby heating the system. The temperature 
after shaking shows resonances as a function of frequency 
[69] . Alternatively, RF spectroscopy can be used [701 [71] . 



VII. CONCLUSIONS AND OUTLOOK 

In this work we analyzed few-body complexes formed 
of dipolar molecules confined in two one-dimensional 
tubes. We studied the stability of these few-body com- 
plexes as a function of dipolar moment, transverse con- 
finement, tilting angle, and statistics of the molecules. 
The complexes extend over both tubes and are bound by 
intertube attraction arising from the long-ranged dipo- 
lar interaction. The extension of the complexes allows 
for an optical, non-destructible detection scheme, which 
enables in-situ probing of the complexes. We show that 
the intensity of the scattered light jumps whenever a new 
complex becomes stable but stays constant otherwise. 

Our work can be extended in various ways. As dis- 



cussed in Sec. IIIF we expect similar results for bi- 
layer systems instead of bitubes. Another extension con- 
cerns possible few-body states for a larger number of 
tubes/layers. In the case of more than two adjacent tubes 
or layers, more complex few-body states become possi- 
ble. In particular, chains of molecules can form with one 
molecule per tube/layer [11 . By tilting the direction 
of the dipoles more than one molecule per layer can be 
bound which can lead to a bifurcation of the tubes. 

In Sec. [IV] we addressed the impact of our results on 
the many-body systems with a finite density in each tube. 
Particularly interesting are fermionic systems. There the 
analog of the dimer physics in the many-body context is 
the BEC-BCS crossover at finite density [33]. In the BEC 
regime the effective entities of the systems are bosonic 
and correspond to tightly bound dimers, while in the 
BCS regime the attraction leads to Cooper pairing but 
still the Fermi-energy is larger than the binding energy 
of the dimer state. However, such a crossover will be 
modified by the presence of trimers. 



We found for the few-body states an extended regime, 
where trimer are stable but dimer do not yet bind 
(forming the 2-2 complex). Translated to the many- 
body regime, this suggests that few-body complexes like 
trimers could be relevant for an intermediate interaction 
regime before the system undergoes Wigner crystalliza- 
tion in the strong interaction limit [23 . These examples 
of many-body physics are relevant for Fermi gases close 
to degeneracy, where the temperature is of the order of or 
smaller than the Fermi energy. However, we note that for 
strong dipole interaction stable complexes will be present 
even in experiments in thermal gases away from degener- 
acy. In particular the proposed detection scheme should 
also work on thermal clouds. 
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Appendix A: Relative coordinates and space grid 

The center of mass motion is unaffected by interactions 
and can be separated by introducing relative coordinates. 



H = H CM (u ) + H Rel ({u}) (Al) 

H Rel ({u}) = ^ (E^ 1 dl 3 + U V({u})) , 

where {u} = (ui, Wjv-i) collectively denotes the TV — 1 
relative coordinates. There is freedom in choosing rela- 
tive coordinates and we choose a transformation matrix, 
M, 



/ u ) 




( x \ 




\-M\ 




\ WJV-1 ) 




V %N-1 I 



such that the kinetic energy is proportional to the unit 
matrix. We chose the following forms of M relating {u} 
with the original coordinates for the different few-body 
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complexes 



M D = l/y/2 



1 1 
-1 1 



l/y/3 l/y/3 l/y/3* 
M Tr = ( -2/y/E l/x/6 l/\/6 
-1/^2 l>/2 



2-2 



/ 1/2 1/2 1/2 1/2 \ 

1/2 1/2 -1/2 -1/2 

-l/y/2 l/y/2 

V -I/a/2 ly/2 J 



(A2) 



(A3) 



(A4) 



M 



1-3 



V 



1/2 1/2 1/2 1/2 \ 

-a/3/2 1/2^3 l/2>/3 1/2^3 

V^/V^ -1/V6 1/a/6 

-I/a/2 1/a/2 J 



(A5) 



The wave functions have the following symmetries: 

^Tr(u U U 2 ) = Slp Tr {u 1 ,-U 2 ) 

^2-2(^1,^2,^3) = si/j F i(u 1 ,-u 2 ,us) 

= S1p F1 (u!,U2, -Us) 
= ^2-2(^1,^3,^2) 
^1-3(^1,^2,^3) = Sl/> F2 (ui,U2, -U 3 ) 
-u 2 + V3u 3 V3U2+U3 s 



2—VSu 3 



-V3u 2 +u 3 
2 



Corresponding to the symmetries we use a square and 
cubic lattice for trimer and 2-2, while for 1-3 we choose a 
hexagonal lattice in the u 2 — us plane and a regular grid 
in ui plane. The 2-2 ground state will have the symmetry 

p = i. 

The Laplace operator in square lattices with lattice 
constants a\ and a 2 takes the form 



A/(x) = J2 I/( x + a &) ~ 2 /( x ) + /( x " Oiei)] l a 1 

i=l 

(A6) 



with ei = (1,0) and e 2 = (0, 1). In the hexagonal lattice 
with lattice constant a it is given by 



A/(x) = - V [/(x + oei) - 2/(x) + /(x - a^)] /a 2 
J i=i 

(A7) 



with ei = (—1,0), where e 2 = ^(\/3, — 1) and e2 = 
|(V3,1). 
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